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In the Arctic, earlier snowmelt and longer growing seasons due to warming have been hypothesized to increase 
vegetation productivity. Using the Normalized Difference Vegetation Index (NDVI) from both field and satellite 
measurements as an indicator of vegetation phenology and productivity, we monitored spatial and temporal pat- 
terns of vegetation growth for a coastal wet sedge tundra site near Barrow, Alaska over three growing seasons 
(2000-2002 ). Contrary to expectation, earlier snowmelt did not lead to increased productivity. Instead, productiv- 
ity was associated primarily with precipitation and soil moisture, and secondarily with growing degree days, 
which, during this period, led to reduced growth in years with earlier snowmelt. Additional moisture effects on 
productivity and species distribution, operating over a longer time scale, were evident in spatial NDVI patterns 
associated with microtopography. Lower, wetter regions dominated by graminoids were more productive than 
higher, drier locations having a higher percentage of lichens and mosses, despite the earlier snowmelt at the 
more elevated sites. These results call into question the oft-stated hypothesis that earlier arctic growing seasons 
will lead to greater vegetation productivity. Rather, they agree with an emerging body of evidence from recent 
field studies indicating that early-season, local environmental conditions, notably moisture and temperature, 
are primary factors determining arctic vegetation productivity. For this coastal arctic site, early growing season 
conditions are strongly influenced by microtopography, hydrology, and regional sea ice dynamics, and may not 
be easily predicted from snowmelt date or seasonal average air temperatures alone. Our comparison of field to 
satellite NDVI also highlights the value of in-situ monitoring of actual vegetation responses using field optical 
sampling to obtain detailed information on surface conditions not possible from satellite observations alone. 

© 2012 Elsevier Inc. All rights reserved. 


1. Introduction 

Global circulation models predict that current atmospheric forcing 
will lead to significant warming of the terrestrial surface in the next 
few decades, with the most rapid warming expected in northern 
latitudes (ACIA, 2005; IPCC, 2007). Long-term temperature trends 
confirm this geographic pattern of warming: the North Slope of Alaska 
has warmed rapidly, with an increase of approximately 2-3 °C reported 
for northern Alaska, over recent decades (Chapin et al„ 2005; Hinzman 
et al., 2005; Overland et al„ 2002). One possible consequence of north- 
latitude warming is a longer growing season, coupled with enhanced 
greening and productivity of northern latitudes. Indeed, satellite 
observations (Bhatt et al„ 2010; Goetz et al„ 2005; Myneni et al., 
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1997) suggest that earlier growing season and enhanced greening are 
already occurring in many arctic regions, although independent field 
confirmation of satellite observations has been largely lacking, with 
the exception of Bhatt et al. (2010). Historical records of vegetation 
in Alaska's North Slope provide additional evidence of enhanced 
north-latitude greening in the form of increased vegetation growth 
(Bhatt et al„ 2010; Chapin et al„ 2005; Sturm et al„ 2001). Because 
warming has been particularly rapid in Alaska (Chapin et al., 2005; 
Hinzman et al„ 2005), regime changes in this region may be indicative 
of how other northern terrestrial ecosystems will change in response 
to similar forcing mechanisms. 

Of particular interest and concern is the response of arctic landscapes 
to altered climate and the potential for shifting biosphere-climate feed- 
backs. It is generally assumed that warming of northern latitudes will 
enhance photosynthetic productivity (hypothesis #1), and the reported 
increased greening and vegetation cover in northern latitudes (Bhatt 
et al., 2010; Goetz et al„ 2005; Myneni et al„ 1997; Sturm et al„ 2001) 
provides evidence for enhanced productivity. However, experimental 
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studies and direct measurements of gas exchange or carbon balance in 
recent decades suggest that tundra ecosystems can be net sources of 
carbon to the atmosphere, due to enhanced ecosystem respiration with 
warming (hypothesis #2, Oechel et al., 1993; Mack et al., 2004; Schuur 
et al., 2008). Clearly, both photosynthetic and respiratory processes can 
be stimulated by warming (Oberbauer et al., 2007), and the balance of 
these processes determines the net carbon flux (Huemmrich et al., 
2010a, b; Jahn et al., 2010). Despite many local observations of enhanced 
photosynthesis and respiration with increasing temperature, we know 
little about how the net biosphere-atmosphere carbon balance may be 
changing for the Arctic as a whole, primarily because of the vast areas 
and different vegetation types involved, most of which are relatively 
inaccessible. This sampling challenge suggests that remote sensing can 
play a key role in resolving the carbon balance. However, to be effective, 
retrievals from remote sensing must be validated using field observa- 
tions that are well-grounded in a solid ecophysiological understanding 
of the basic ecosystem processes affecting the carbon balance. We 
demonstrate that field optical sampling (proximal remote sensing) can 
play an important role in this validation. 

Multi-factor studies of climate change effects on terrestrial ecosys- 
tems demonstrate that carbon balance cannot always be predicted by 
changes in temperature alone, but is determined by a complex interplay 
of many factors, all of which may be changing simultaneously (Field et al„ 
1992; Mooney et al., 1991). Likely factors include moisture availability, 
nutrient availability, species composition, and atmospheric carbon diox- 
ide levels, all of which interact with temperature. The productivity and 
carbon balance of arctic vegetation are known to be particularly respon- 
sive to moisture (Bliss et al., 1984; Chapin et al., 1995; Oechel et al., 1993) 
and nutrients (Shaver et al., 2001 ), both of which can vary over short dis- 
tances and time scales and are likely to change with climate. Small-scale 
variations in factors that determine productivity and composition make it 
difficult to use coarse-scale satellite data to understand or quantify arctic 
vegetation changes in response to regional climate variation. 

Temperature and moisture are key variables determining the 
carbon balance of arctic tundra (Huemmrich et al., 2010a, b; Jahn et al., 
2010). Tundra ecosystems are particularly sensitive to prevailing mois- 
ture conditions, because production is often moisture-limited in tundra 
ecosystems (Bliss et al., 1984; Huemmrich et al., 2010 a,b). 

For arctic tundra, photosynthesis, respiration, and production tend 
to increase with moisture, but these patterns can reverse when the 
surface becomes fully inundated (Dagg & Lafleur, 2011). Both respira- 
tion and photosynthesis are also strongly temperature-dependent, but 
exhibit different temperature optima that vary with species and growth 
conditions (Berry & Bjorkman, 1981; Field et al„ 1992). In northern 
ecosystems, temperature and moisture are often linked. For example, 
recent warming has been accompanied by increasing summer thaw 
depth, which often results in localized draining and drying of surface 
soils (Jorgenson et al., 2006; Smith et al., 2005). Further complicating 
the prediction of carbon balance is the fact that ecosystem respira- 
tion is composed not only of the respiration of recently-fixed carbon 
(i.e. autotrophic respiration from vegetation), but also involves 
heterotrophic respiration of ancient carbon deposits in deeper 
tundra layers (Jahn et al„ 2010; Schuur et al., 2008). 

While many studies suggest that northern warming will enhance 
vegetation productivity and carbon uptake by arctic ecosystems 
(e.g. Goetz et al., 2005; Groendahl et al., 2007; Myneni et al., 
1997), an alternate consideration is that earlier snowmelt in 
northern latitudes associated with a warming climate may lead to 
drier growing seasons, resulting in reduced photosynthetic carbon 
uptake. In recent years, large regions of the Arctic have undergone 
surface hydrological changes resulting from a warmer, deeper active 
layer (Smith et al., 2005). Tundra in the North Slope of Alaska may 
be undergoing a regime change towards drier conditions (Mack 
et al., 2011). In Barrow, recent studies have reported vegetation 
cover changes consistent with surface drying and reduced summer 
net primary production (Lara et al., 2012; Lin et al„ 2012). These 


recent observations, combined with a history of reports documenting 
drought limitations to photosynthesis and productivity in arctic vegeta- 
tion (Bliss et al., 1984; Huemmrich et al., 2010a&b; Oechel et al., 1993), 
suggest that drought restraints could override any warming benefits to 
vegetation growth. Mosses, which account for a large fraction of arctic 
vegetation cover, are particularly sensitive to such drying and rapidly 
shut down photosynthesis during dry periods (McNeil & Waddington, 
2003; Silvola, 1990). Thus, it is particularly likely that tundra productiv- 
ity and carbon balance would be strongly influenced by the particular 
combination of temperature and moisture experienced during the 
short growing season. 

A challenge to understanding northern latitude ecosystem responses 
to rapid climate change has been the lack of sustained observations of 
ecosystem processes. There is a strong need for consistent, repeatable 
sampling of ecosystem properties and processes across the Arctic. The 
satellite record offers some promise, but also presents several technical 
difficulties that are particularly profound in northern latitude sampling. 
These challenges include persistent cloud cover, low sun angles, errors 
arising from temporal compositing methods, the lack of standard band 
definitions, the poor resolution of fine-scale surface features due to 
large pixel sizes, and the confounding effects of water (variably present 
as liquid water, snow, or ice, each with its own characteristic optical 
properties) (Stow et al., 2004). All of these factors complicate the charac- 
terization of surface properties, and can lead to ambiguous conclusions 
regarding changes in vegetation cover and physiological activity for 
arctic regions. For example, satellite detection of apparent surface 
“greening" attributed to vegetation growth can be caused by a reduction 
in snow cover (associated with an earlier snowmelt) or by a reduction in 
the area of standing water (associated with increasing thaw depth). 
Given the strong role remote sensing has played, and is likely to play, 
in understanding biotic responses to climate change, and given the 
technical challenges described above, there is a particular need for 
field-based optical measurements of the tundra surface that can inde- 
pendently validate satellite observations. Because weather and climate 
may have different local effects depending upon microclimate, the ability 
to distinguish microtopographic effects on surface conditions at a high 
temporal resolution may be particularly useful in revealing underlying 
functional processes contributing to ecosystem change. 

To meet the need for repeatable, accurate surface measurements 
that characterize biospheric responses, we implemented a field 
sampling protocol for surface reflectance involving repeated transect 
sampling with a field spectrometer. The fine spatial, temporal, and 
spectral scales of this sampling allowed detailed observations of the 
response of tundra productivity to seasonal and interannual changes 
in weather conditions as a proxy for changing climate. By sampling at 
the surface-atmosphere interface (under the atmosphere rather than 
through the atmosphere), we were able to avoid the confounding 
effects of cloud cover that often leads to ambiguous retrievals of surface 
reflectance. This sampling method was implemented over three succes- 
sive growing seasons to reveal the detailed response of an arctic tundra 
ecosystem to a variety of factors (microtopography, species composi- 
tion, temperature, and moisture). 

The primary goals of this study were to better quantify the spatial 
and temporal dynamics of tundra productivity in response to changing 
temperature and moisture regimes, and to test the hypothesis that an 
earlier growing season leads to greater vegetation productivity. In this 
case, NDVI (an optical index of vegetation greenness), was used as a 
proxy for productivity based on observations of significant correlations 
between NDVI, above-ground biomass, and leaf area index for nearby 
tundra vegetation (Goswami, 2011). 

2. Methods 

All field optical measurements were conducted along a 100-m tran- 
sect (center lat-lon coordinates: 71.32202265 and — 156.60297), located 
approximately 300 m southeast of the National Oceanic and Atmospheric 
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Administration (NOAA) Earth System Research Laboratory (ESRL) Global 
Monitoring Division (GMD) baseline observatory near Barrow, AK (BRW). 
Our research site is included in the SpecNet network (Spectral Network — 
Gamon et al., 2006b, 2010). The dominant vegetation of this area is wet 
sedge tundra, typical of coastal arctic, containing “patterned ground" 
with characteristic troughs and high-centered polygons, a landscape pat- 
tern common near Barrow, Alaska (Brown et al., 1 980). Dominant species, 
determined by sampling meter-square plots along the transect in August 
2001, are listed in Table 1. 

2.1. Meteorological data 

Daily average broadband solar irradiance was measured at the adja- 
cent NOAA GMD (BRW) site. Air temperature at a height of 2 m was also 
measured at BRW and expressed as daily averages. Precipitation data 
were measured at the National Weather Service's (NWS) Barrow sta- 
tion, and were obtained from http://www.climate.gOv/#dataServices, 
and expressed as daily totals of water-equivalent precipitation. To 
help visualize inter-annual differences, temperature was also expressed 
as growing degree days (cumulative values above 0 °C), and precipita- 
tion was also expressed as cumulative totals, using day 131 as a starting 
date. 

2.2. Soil moisture 

Soil moisture was sampled at the “Barrow 1" site (Soil Survey Staff, 
Natural Resources Conseivation Service, United States Department of 
Agriculture, Soil Climate Research Station Data, accessed May 24, 
2012, from http://soils.usda.gov/suivey/smst/). The sampling location 
(71.32238923, — 156.6107969 W) was approximately 250 m from the 
optical sampling transect on a similar landscape and elevation. Soil 
moisture was sampled at 3 depths (5, 20, and 30 cm) with 4 replicates 
during the study period (2000-2002) using 12 dielectric constant 
soil moisture/temperature sensors (Hydra Probe, Stevens Vitel, Inc., 
Chantilly, VA, USA) attached to dataloggers (Campbell Scientific, 
Logan, UT, USA). Data from the 12 sensors were averaged to create a 
single average soil moisture value, expressed as a unitless volumetric 
fraction of water (vfw), with units of m 3 water per m 3 soil. Note that 
since the dielectric constant of water increases as water thaws, and is 
also responsive to liquid soil water content (Hydra Soil Moisture 
Probe User's Manual, 1994, Stevens Vitel, Inc., Chantilly VA, USA), this 


Table 1 

Common species and genera found along the sampling transect by category (lichen, 
moss, and vascular plant). 


Mosses 

Lichens 

Vascular plants 

Aulacomnium turgidum 

Alectoria nigricans 

Graminoids 

Dicranum elongatum 

Cetraria cucullata 

Arctagrostis latifolia 

Drepanocladus revolvens 

Cetraria nivalis 

Carex aquatilis 

Drepanocladus uncinatus 

Cladina rangiferina 

Dupontia fisherii 

Polytrichum sp. 

Cladonia chlorophaea 

Eriophorum angustifolium 

Sphagnum sp. 

Cladonia scabriuscula 

Eriophorum scheuzerii 


Cladonia squamosa 

Eriophorum vaginatum 


Dactylina arctica 

Luzula confusa 


Ochrolechia frigida 

Poa arctica 


Pertusaria sp. 

Dwarf-shrubs 


Sphaerophorus globosus 

Salix pulchra 


Thamnolia subuliformis 

Salix rotundifolia 
Saxifraga spp. 
Vaccinium vitis-idaea 
Cochlearia officinalis 
Forbs 

Potentilla hyparctica 
Ranunculus nivalis 
Senecio atropurpureus 
Stellaria edwardsii 
Petasites frigidus 
Papaver macounii 


method provides an estimate of soil water available to plants during 
the early growing season. 

2.3. Spectral reflectance 

Spectral reflectance measurements were made at each meter 
along a 100-meter transect running east-west across the landscape 
(center lat-lon coordinates were 71.32202265 and —156.60297, 
respectively). In most cases, reflectance measurements employed a 
portable field spectrometer with two detectors (UniSpec DC, PP 
Systems, Amesbury MA), one sampling surface radiance, and the 
other sampling sky irradiance. These two signals were compared 
over a white reference standard (Spectralon, LabSphere, North 
Sutton, New Hampshire, USA), yielding a cross-calibration used to 
calculate reflectance (see Gamon et al., 2006a for details). By simulta- 
neously correcting each surface radiance measurement by the 
prevailing irradiance, this dual-detector spectrometer corrected for 
changing light conditions due to the frequent cloud cover in this 
region. Most days during the sampling period were partly cloudy 
or fully overcast, rendering optical satellite observation difficult 
(see below). When the dual-detector spectrometer was not available, 
reflectance was sampled by a single-detector spectrometer (UniSpec, 
PP Systems, Amesbury MA), and reflectance calculated by dividing the 
surface radiance against the radiance of the same white reference stan- 
dard. For the spectral regions reported here, reflectance measurements 
from the two instruments were statistically indistinguishable on the 
basis of a Student’s t-test. Consequently, single-detector spectrometer 
measurements were substituted for those of the dual-detector spec- 
trometer when the dual-detector was not available. 

To improve sampling precision, and to enable point-by-point 
seasonal comparisons within a given year, most measurements 
during the growing season were conducted by mounting the dual- 
detector spectrometer on a “tram system" consisting of a mobile 
cart on a track mounted on a series of tripods raised above the tundra 
surface (Gamon et al., 2006a; Goswami et al., 2011). In this configura- 
tion, the downward-looking foreoptic was mounted on a boom 
positioned approximately 2.5 m above the tundra surface. At this 
range, each spectral reflectance measurement sampled a surface 
area approximately 0.8 m in diameter. Sampling from the tram was 
not always possible, in part because local safety concerns required 
that the tram tracks be removed for the winter and reinstalled at 
the start of the next season. To gather data prior to snowmelt and 
tram deployment each season, measurements were made in the iden- 
tical tram positions using a manual sampling method that replicated 
the tram sampling positions using permanent wooden markers and 
a 100-meter sampling tape. To do this, we walked the transect taking 
a sample at each meter mark, using a leveled, hand-held mast and 
boom that replicated the tram sampling position. In 2000 and 2001, 
the tram was installed just around the time of spring snowmelt and 
removed at the end of the growing season, thus avoiding snowdrifts 
that would have led to artifacts in the snowmelt date. In 2002, the 
entire sampling was conducted manually because the tram system 
was not available during that year. Transect reflectance measure- 
ments were sampled 2-3 times a week during each growing season. 

Reflectance spectra were interpolated to yield reflectance for 1-nm 
wavelength intervals (“bands") between 400 and 1100 nm. The 
Normalized Difference Vegetation Index (NDVI), a widely used measure 
of arctic vegetation greenness and productivity (Stow et al„ 2004), was 
then calculated from these interpolated spectra using the following 
formula: 

NDVI = (Pnir Pred)/ (Pnir + Pred) (1) 

where p N1R and p RED indicate reflectance in the near-infrared (NIR) and 
red bands, respectively. In this study, we used 670 nm for the red 
band, and 800 for the NIR band. Because NDVI measurements were 


JA. Gamon et al. / Remote Sensing of Environment 129 (2013) 144-153 


147 


nondestructive, they allowed continuous sampling of the same regions 
during the growing season. Note that some reports (e.g. Hope et al., 
1993) have reported variable NDVI-biomass relationships for different 
arctic vegetation types, so we examined this issue for vegetation of 
our site. Calibration data collected from a nearby site confirmed a strong 
correlation between NDVI and above-ground, biomass for several com- 
mon study species (R 2 = 0.83 for dominant sedges and R 2 = 0.78 for all 
species combined), confirming the use of NDVI as a useful proxy of 
above-ground vegetation productivity for this ecosystem, particularly 
when vegetation composition was not markedly changing, as was the 
case during this study. Harvest results and NDVI calibrations have 
been reported elsewhere (Goswami, 2011) so they are not presented 
in further detail here. 

NDVI was calculated for each meter mark along the transect. In 
some cases, these 100 samples for each date were averaged to create 
a single mean value for the transect on a given sampling date. These 
transect-averaged values, also termed “big-pixel" values for their 
potential to represent a single pixel for a low-resolution satellite 
sensor, were used to more clearly depict seasonal and year-to-year 
trends for the site as a whole. To examine spatial and temporal trends 
at a finer resolution and enable tests of statistical significance, we also 
calculated NDVI data from each of the (approx. 0.8 m diameter) 100 
sampling locations, termed “small-pixel values.” 

We also measured reflectance at 51 plots, each about 15 cm in 
diameter, representing nearly pure areas of representative tundra 
species. These plots were located within 15 m of the transect and 
were measured multiple times in 2001 to describe seasonal changes 
for representative species from each major cover type (Table 1 ). 

2.4. MODIS NDVI data 

NASA MODIS TERRA NDVI data (MOD13Q1, subsetted land prod- 
ucts, collection 5) were obtained from the Oak Ridge National Labora- 
tory Distributed Active Archive Center (ORNL DAAC, Oak Ridge, 
Tennessee, USA) available at http://daac.ornl.gov/MODIS/modis.html 
(accessed July 8, 2012). NDVI data were downloaded as a single, 
250-m pixel matching the location of the optical sampling transect 
(71.32202265, —156.60297), and data were sorted by quality flag 
values into “good data” (reliability 0) and other quality categories 
(reliability 1-3) for comparison with ground NDVI. Single-pixel 
NDVI values were selected to minimize possible contamination from 
standing water, roads, and buildings in nearby pixels. 

2.5. Snowmelt dates 

Snowmelt was determined from surface reflectance and is 
reported two ways. NDVI values were very sensitive to snow cover, 
and were compared to field observations of snow cover to evaluate 
the start and stop of snowmelt along the transect. The start of 
snowmelt, when some bare ground became visible, was defined as 
the date when the NDVI first became positive. The end of snowmelt 
(i.e. snow-free surface, based on field observations) was defined as 
the date when NDVI reached 0.3. We also used the NOAA GMD 
BRW snowmelt dates, measured approximately 300 m from the 
tram transect. BRW snowmelt was determined as the day of year 
when the average broadband albedo dropped below a threshold 
value of 0.3 (Stone et al., 2002). 

3. Results 

3.1. Meteorological patterns 

Downwelling solar irradiance followed a similar cycle in each of 
the three years, increasing towards summer solstice, then declining 
(Fig. 1A). Clearly noticeable is a marked increase in day-to-day varia- 
tion during the summer growing season (roughly days 150-270) due 


to increased cloudiness following snowmelt. Cloud cover at Barrow 
tends to increase during the transition from spring to summer as 
the snow and sea ice in the adjacent Beaufort Sea melts and releases 
moisture into the atmosphere. Arctic stratus clouds form and are 
most common during late summer, when sea ice is at its minimum 
(Stone et al., 1996). Despite considerable intra-annual variation due 
to cloud cover (Stone, 1997), differences in cumulative irradiance 
(Fig. IB) or average irradiance (Table 2) during the greenup period 
(DOY 131-220) between years were not significant. 

Air temperature trends followed similar patterns in 2000 and 
2001, increasing steadily towards day 160 (June) then hovering 
between 0 and 10 °C for most of the summer (Fig. 2A). In 2002, air tem- 
perature was unusually warm early in the season (days 130-150), lead- 
ing to an early snowmelt (days 141-157 along the transect; Table 2) 
then was slightly cooler for most of the subsequent months (days 
160-190) than the previous two years (Fig. 2A). Transect snowmelt 
dates fell within the range of values obtained from the nearby BRW 
station (300 m away), and had the same rank order, with each year 
having an earlier melt than the previous year (Table 2). 

These temperature patterns can best be seen when plotted as 
“growing degree days" (summed daily mean values above 0 °C) 
(Fig. 2B). Interestingly, after an initial rapid spring increase, growing 
degree days for 2002 lagged behind that of the previous two years 
for much of the growing season, and this delay was particularly evi- 
dent during days 170-195 (Fig. 2B). From days 195 to 240, covering 
the peak of the growing season, growing degree days were similar 
in 2001 and 2002, and both lagged the values for 2000 (Fig. 2B). 

Growing season precipitation was generally low, and the pattern 
and amount of precipitation varied among years, with the greenup 
period of 2000 having the highest precipitation, and that of 2002 
the lowest (Table 2). In 2002, precipitation was low during the first 
half of the growing season (particularly between days 180 to 220, 
Fig. 3A). The lower early-season precipitation in 2002 was more 



Fig. 1. Daily average downwelling solar irradiance for 2000-2002 for the entire year 
(top) and cumulative irradiance for the active growing season (bottom). Symbols in 
panel B indicate NOAA snowmelt dates (see Table 2) for each year. Peak season was 
determined by maximum NDVI (approximately day 220). 
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Table 2 

Snowmelt date, meteorological variables, and NDV1 in 2000-2002. Snowmelt dates were 
determined from the study transect and from the nearby NOAA GMD BRW station (see 
Methods). Average daily irradiance, average daily temperature, average daily precipitation, 
total precipitation, and volumetric soil water content (vfw, expressed as m 3 water per m 3 
soil) are shown during the greenup phase (DOY 131-220). Growing degree days at peak 
season (day 220), peak NDVI, and seasonally-integrated NDVI are shown for 2000, 2001, 
and 2002. Significant differences (p<0.05) in average meteorological values and peak 
NDVI values between years are indicated by different superscripts. Seasonally-integrated 
NDVI was calculated from snowmelt to Julian day 255, using the snow-free NDVI value 
for this site (0.3) as a baseline. 



2000 

2001 

2002 

Transect snowmelt period (day of year) 

163-168 

162-165 

141-157 

NOAA snowmelt date (day of year) 

165 

162 

144 

Ave. daily irradiance (Wm -2 ) 

214 a 

22T 

221 a 

Ave. daily temperature (°C) 

0.41 ab 

-0.69 b 

0.76 a 

Growing degree days at peak season 

194 

142 

148 

Ave. daily precipitation (mm) 

1.06 a 

0.77 a 

0.44 b 

Total precipitation (mm) 

95 

70 

39 

Ave. volumetric soil water content (vfw) 

0.1 89 a 

0.1 84 a 

0.179 3 

Peak NDVI 

0.585 a 

0.573 b 

0.526 c 

Integrated NDVI 

16.7 

13.4 

12.2 


clearly evident when expressed as cumulative precipitation (Fig. 3B), 
demonstrating a clear lag by peak season in 2002 relative to the 
other years. 

Following snowmelt, soil moisture closely followed precipitation 
patterns, particularly between days 180 and 220 leading up to the 
peak season. In 2000, a large rainfall event on day 187 followed by 
several smaller events (Fig. 3 A and B) caused a steady increase in 
both cumulative precipitation and soil moisture towards the peak 
season (Fig. 3B and C). In 2001, the lack of rain between days 180 
and 200, followed by several rain events between days 200 and 220, 
caused a slight decline in soil moisture (days 180-200), followed by 
an increase to 2000 levels by mid-season (days 200-220, Fig. 3C). In 




Fig. 2. Daily average temperature (top) or growing degree days (cumulative values 
above 0 °C, bottom) for 2000-2002. Symbols in panel B indicate NOAA snowmelt 
dates (see Table 2) for each year. Peak season was determined by maximum NDVI 
(approximately day 220). 


2002, after some initial rain between days 170 and 181, very little 
additional rain fell in the first half of the growing season, causing 
both cumulative precipitation and soil moisture to lag behind the 
values for the previous two years (Fig. 3B and C). 


3.2. Microtopographic effects on surface reflectance 

Spatial and temporal patterns of surface reflectance along the tran- 
sect revealed clear microtopographic effects on the timing of snowmelt 
and vegetation growth. A spatial-temporal plot of NDVI from 2001 is 
shown in Fig. 4. Areas to first emerge from snow (e.g. meters 18 and 
65) were the local high spots in the landscape. By the peak of the grow- 
ing season (late July to early August) these high spots that first became 
snow-free had relatively low NDVI values. By contrast, low-lying 
regions slow to emerge from snow reached the highest NDVI values 
by mid-season (e.g. meters 10 and 82). 





Fig. 3. Daily precipitation totals (A), cumulative precipitation (B), and column- 
averaged soil water content expressed as volumetric fraction of water (C) for the 
2000-2002 growing seasons. Symbols in panels B and C indicate NOAA snowmelt 
dates (see Table 2) for each year. Peak season was determined by maximum NDVI 
(approximately day 220). 
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These microtopographic effects can be seen more clearly in a “delta 
NDVI" plot that shows the seasonal increase in NDVI from snowmelt 
to mid-season (Fig. 5). By referencing mid-season NDVI to the early- 
season, snowmelt value (0.3), this figure removes NDVI artifacts 
due to snowmelt, and highlights the NDVI changes due to vegetation 
growth only. The total difference in elevation along the 100-m transect 
was small (1 m), with local variations typically much less than this 
(dashed line, Fig. 5, bottom panel), yet these minor elevation differ- 
ences had a profound effect on species composition and productivity 
as indicated by the seasonal increase in NDVI. Slight local depressions 
(typically 10-20 cm in depth) tended to have higher moisture content, 
often containing standing water for short periods. These wet locations 
were dominated by vascular plants, particularly graminoids (Fig. 5, 
top panel). These were also the regions with the highest seasonal pro- 
ductivity gains, measured by the mid-season NDVI increase (Fig. 5, bot- 
tom panel). By contrast, higher elevation microsites had a greater cover 
of mosses and lichens, were generally drier, and showed less seasonal 
change in NDVI (Fig. 5). 

3.3. Seasonal and interannual patterns of surface reflectance 

Comparison of “big-pixel" NDVI values over three successive 
growing seasons (2000-2002) revealed inter-annual differences in 
snowmelt, growing season length, and productivity, with each year 
having an earlier start to the growing season, yet reaching a lower 
maximum NDVI (Fig. 6). In 2000, complete snowmelt (defined as 
the date when average NDVI reached the snow-free value of 0.3) 
occurred on day 168, and NDVI peaked at a value of 0.585 (day 217). 
In 2001, complete snowmelt occurred three days earlier on day 165, 
yet the production lagged behind the previous year, with NDVI reaching 
a maximum value of 0.573 (day 216). This peak value was significantly 
different from the peak value in 2000 (p = 0.002, paired samples t-test). 
In 2002, complete snowmelt along the transect occurred on day 157, 
eleven days earlier than in 2000. At the nearby NOAA site, 2002 was 
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Fig. 5. Vegetation cover (top panel), and seasonal increase in NDVI from snowmelt to 
mid-season in 2000 (ANDVI) and local elevation (bottom panel) plotted for each 
meter along the tram transect. Vegetation cover was sampled in early August 2001 
and represents the cover at peak-season. Total ground cover exceeds 100% due to 
multiple canopy layers (e.g. vascular plants overtopping moss or lichens). ANDVI 
was determined as the difference between day 220 (maximum NDVI) and day 166 
(immediately after snowmelt). 
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Fig. 4. Spatial and temporal patterns of “small-pixel" NDVI measured along the tram 
transect in 2001 for the first part of the growing season. Snow-covered values having 
low NDVI are shown in blue. The rapid rise in NDVI between days 162 and 165 is caused 
by snowmelt (Table 2), and the more gradual rise following snowmelt is due to vegetation 
growth. (For interpretation of the references to color in this figure legend, the reader is 
referred to the web version of this article.) 


determined to be the earliest snowmelt on record to date (snowmelt 
records for this site date to the early 1940s, Stone et al„ 2005; see also 
http://www.esrl.noaa.gov/gmd/grad/snomelt.html). Following this re- 
markably early snowmelt, production lagged behind that measured in 
the previous two years, with NDVI reaching a maximum NDVI value 
of only 0.526 (day 222). Despite this unusually early start to the grow- 
ing season, the peak NDVI value in 2002 was significantly lower than 



Fig. 6. Seasonal course of the Normalized Difference Vegetation Index (NDVI), an index 
of “greenness” and plant productivity, during the 2000, 2001, and 2002 growing 
seasons (samples collected 2-3 times per week). In this case, each “big pixel" NDVI 
value was calculated from reflectance averaged from the entire transect (100 sampling 
points). Arrows indicate dates of complete snowmelt along the transect based on 
visual observation, which coincided with NDVI values exceeding 0.3. Symbols indicate 
NOAA snowmelt dates, all of which fell within the transect snowmelt period (NDVI 
between 0 and 0.3). See Table 2 for a summary of snowmelt dates. 
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the peak NDVI values in either of the two previous years (p<0.001, 
paired samples t-test, Table 2). Integrated NDVI values followed a 
similar trend, with the lowest integrated NDVI values occurring in 
2002 (Table 2). 

Reflectance for the three major cover types - mosses, lichens, and 
vascular plants - revealed contrasting NDVI patterns over the 2001 
growing season (Fig. 7). The NDVI for vascular plants closely followed 
the NDVI for the transect as a whole — increasing towards mid-season 
(approximately day 220), then declining. On the other hand, moss NDVI 
initially declined, then gradually increased. This increase followed a 
period of rainfall starting on day 199 in 2001 (Fig. 3). Lichen NDVI 
followed a similar decline and increase after the rainfall events. These 
results reflect the strong sensitivity of physiological activity to moisture, 
and indicate that the growing-season NDVI patterns at this site are 
primarily driven by the growth and senescence of the vascular plants. 

3.4. Comparison of ground to satellite NDVI 

MODIS 16-day composite NDVI values were compared to ground- 
measured NDVI values from Fig. 6, revealing several key points. In 
all three years, MODIS NDVI values showed a similar seasonality to 
ground-based values (Fig. 8). The highest quality MODIS values 
(solid circles, Fig. 8) tended to be slightly higher than the correspond- 
ing ground NDVI values. However, MODIS composite NDVI could not 
accurately resolve snowmelt, sometimes appearing to lead (Fig. 8A) 
and other times to lag (Fig. 8B and C) actual snowmelt. Similarly, 
MODIS was unable to pinpoint the date or maximum value of 
peak NDVI, sometimes yielding maximum values that were later 
than (Fig. 8A and B) or earlier than (Fig. 8C) the maximum values 
measured on the ground. 

4. Discussion 

Our findings contradict the hypothesis predicting greater produc- 
tion under longer growing seasons. We expected earlier spring melt 
and a longer snow-free growing season to result in greater vegetation 
growth, as had been suggested for other northern latitude ecosystems 
(e.g. Bhatt et al., 2010; Goetz et al., 2005; Myneni et al., 1997). Instead, 
our measurements indicated an unexpected and progressive decline in 
vegetation growth (sampled as NDVI) as snowmelt advanced each 
year (Fig. 6), with the lowest midseason NDVI attained in 2002, the 
year with the earliest snowmelt on record (Stone et al., 2005). Clearly, 
factors other than the timing of snowmelt determined the overall 
seasonal productivity during these years at this site. 

Reduced production with earlier snowmelt was also apparent in 
the microtopographic analysis, where sites that emerged earlier 



Day of year 


Fig. 7. Seasonal progression of NDVI for three dominant vegetation cover types (vascular 
plants, mosses and lichens) along the transect in 2001 . Also shown for comparison is the 
average NDVI pattern for the entire transect from Fig. 6. 




Fig. 8. Comparison of MODIS composite NDVI (circles and dashed lines) to 
ground-measured NDVI (solid lines) for 2000 (A), 2001 (B) and 2002 (C). Solid circles 
indicate MODIS NDVI data having the highest quality flag (reliability = 0); open circles 
indicate points with lower quality flags ( reliability = 1, 2 or 3). MODIS values were 
plotted according to actual measurement dates (not nominal composite dates) to facil- 
itate comparison with ground measurements. 

were higher and drier, and tended to have less productivity, mea- 
sured as lower mid-season NDVI values (Figs. 4 and 5). While the 
mechanisms and time scales of the two responses are different, both 
are driven by moisture. Reduced NDVI in successive years was caused 
by reduced vegetation growth due to moisture limitations. On the 
other hand, the spatial differences in NDVI due to microtopographic 
variation represent the result of a long-term vegetation response to 
moisture patterns. These observations, while contrary to simple pre- 
dictions of enhanced productivity with earlier snowmelt, are actually 
consistent with ecophysiological studies linking tundra production to 
moisture and temperature conditions. Our findings most likely reflect 
water limitations to production at this site, and are consistent with a 
body of literature that suggests that the productivity of arctic tundra 
is often moisture-limited (Bliss et al., 1984; Chapin et al., 1995; 
Huemmrich et al., 2010a, b; Oechel et al., 1993). Other arctic field 
studies have reported similar findings, suggesting that arctic vegeta- 
tion productivity is not necessarily correlated with the timing of 
snowmelt, but with environmental conditions prevailing early in 
the growing season prior to peak production (Humpreys & LaFleur, 
2011 ). 

Our observations suggest that conditions during the early growing 
phase, particularly precipitation and soil moisture may be the deter- 
mining factors, rather than irradiance, temperature or growing 
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season length, per se. Of these meteorological variables, precipitation 
and soil moisture best matched the NDVI patterns (Table 2). Despite 
earlier snowmelt in each year, productivity (measured as peak 
NDVI) declined, and this decline closely matched the precipitation 
and soil moisture patterns during the first half of the growing season 
(days 131-220). By contrast, NDVI showed no clear association with 
either solar irradiance or temperature over the same period. Following 
an unusual early-season warming in 2002, cool temperature anomalies 
just after snowmelt may have contributed to the particularly low pro- 
ductivity in that year. The combination of drier and cooler conditions 
in 2002, relative to the previous two years, probably explains the 
lower NDVI for this year. In this study, water limitation appeared to 
override the effects of growing season length and warmer initial 
growing-season temperatures on overall seasonal productivity. 

The importance of moisture limitations to productivity can also be 
seen in the clear relationship between microtopography and NDVI 
for this site, with locally low and wet areas exhibiting the highest 
mid-season NDVI peaks. These peaks were clearly linked to vegetation 
cover, with moisture-loving graminoids (e.g. Carex aquatilis) dominat- 
ing in these locations. By contrast, the elevated microsites, which were 
less than 1 m higher, had a greater percentage of drought-tolerant 
cover types, including prostrate shrubs, mosses and lichens. Our find- 
ings are consistent with other recent reports from this area showing 
contrasting photosynthetic light-use efficiency and productivity asso- 
ciated with microtopography and moisture (Huemmrich et al., 2010, 
submitted for publication). Unlike the interannual differences in 
NDVI that were associated with variation in precipitation and soil 
moisture, these microtopographic effects represent long-term vegeta- 
tion community responses to moisture occurring over much longer 
time spans than the three-year period of this study. Our observations 
of different vegetation cover types with microsite are consistent with 
a large body of literature illustrating arctic vegetation community dif- 
ferences with microtopography (Bliss & Gold, 1994; Brown et al., 
1980; Engstrom et al., 2005; Gamon et al., 2012; Webber, 1978). 

From these observations, we speculate that regional changes in 
hydrology reported for many arctic regions (Jorgenson et al., 2006; 
Smith et al., 2005) could lead to striking changes in surface vegetation, 
productivity, albedo, and other ecosystem properties, all of which could 
exert further feedback effects on the climate and atmosphere through 
changes in surface albedo and carbon balance. These effects would 
include altered productivity and gas exchange in the short term 
(Huemmrich et al., 2010b), and altered vegetation patterns in the 
longer term (Lin et al., 2012; Webber, 1978). Recent studies comparing 
surface conditions in 2010 to those documented in 1972 provide evi- 
dence that some vegetation cover types near Barrow have undergone 
surface drying and declines in net primary production over that time 
period (Lara et al., 2012; Lin et al., 2012). Similarly, regional warming 
and drying recently observed for large regions of the Alaskan North 
Slope (Mack et al., 2011) could be shifting tundra carbon balance 
towards carbon release, despite longer growing seasons. These changes, 
if widespread, could provide a significant positive feedback to a 
warming climate, and may accelerate and enhance regime changes for 
Alaska's Slope tundra ecosystems. 

Our findings contrast with observations from eddy covariance 
measurements made at a nearby site that concluded that 2002 was 
a more productive year than average (Kwon et al., 2006). However, 
eddy covariance measures net carbon exchange, not gross carbon 
exchange, and cannot always be partitioned into component photo- 
synthetic and respiratory fluxes. Kwon et al. (2006) did not report 
gross ecosystem carbon uptake, which would have enabled a more 
direct comparison with our results. In this coastal tundra, ecosystem 
respiration can be a large and variable part of the net carbon flux, 
and the overall balance of carbon uptake and loss is particularly sen- 
sitive to soil moisture (Huemmrich et al., 2010b), which undoubtedly 
varied between years during this study period, leaving it unclear 
what drove the year-to-year variation in the net carbon flux reported 


in Kwon et al. (2006). This inability to partition net fluxes into uptake 
and loss illustrates the difficulty of understanding the driving mecha- 
nisms of flux measurements using eddy covariance methods alone, 
particularly for harsh arctic environments where productivity is low 
and extreme conditions can easily cause instrument failures leading 
to data gaps. This also demonstrates the value of independent valida- 
tion of productivity estimates, and the NDVI sampling method in this 
study provides one way to do this. 

We readily acknowledge that a three-year study is not sufficient to 
establish definitively the relationships between snowmelt date, 
temperature, moisture, and seasonal productivity. Furthermore, the lo- 
cation of this particular site near the coast, where weather is strongly 
influenced by sea-ice patterns, may not be representative of arctic 
sites in general. However, the moisture limitations to productivity 
found here corroborate a well-established history of field and experi- 
mental observations made over many arctic sites, so the general princi- 
ples observed here are likely to be widely applicable. Furthermore, 
vast regions of coastal arctic tundra display similar microtopographic 
features that affect vegetation distribution and related ecosystem 
properties and processes (Brown et al., 1980; Gamon et al., 2012), so 
it is likely that the drought limitations and productivity patterns seen 
in this particular landscape are of larger significance. 

This study was notable in that it included years (2000 and 2001) 
having “typical" snowmelt date, contrasting with another year (2002) 
having the earliest snowmelt on record for the period 1941-2011 
(Stone et al., 2005; see also http://www.esrl.noaa.gov/gmd/grad/ 
snomelthtml). The remarkably early snowmelt in 2002 appeared to 
be associated with regional anomalies in both sea ice and atmospheric 
circulation as described in Stone et al. (2005), supporting the observa- 
tion that the regional synoptic conditions and biological responses of 
Alaska’s North Slope are changing, possibly as a result of global 
warming that has also affected sea ice distributions over the past 
three decades. 

Our observations of a strong rise in NDVI at the time of snowmelt 
(Figs. 4, 6, and 8) illustrate how earlier snowmelt can masquerade as 
greater vegetation production when viewed from a satellite. Both 
snow and standing water are common at this arctic site, and have a 
strong tendency to lower NDVI values relative to vegetated surfaces. 
Consequently, decreases in snow or surface water cover can, in them- 
selves, appear as an apparent greening (visible as a rapid rise to a 
snow-free NDVI value of 0.3 in Figs. 4 and 6). Additionally, the heavy 
cloud cover of many arctic regions during the growing season (Hope 
et al., 1999; Narasimhan Si Stow, 2010; see also Fig. 1), makes recovery 
of accurate surface NDVI values difficult from optical satellites, particu- 
larly when using composited data, as demonstrated by our comparison 
of MODIS NDVI with ground NDVI (Fig. 8). Clearly, composited satellite 
data cannot accurately depict the rapid NDVI transitions during snow- 
melt and early-season greening, and this limitation is further confound- 
ed by cloud cover. In some cases, use of daily MODIS data (rather than 
maximum value compositing) can improve the ability to resolve 
early-season NDVI transitions (Narasimhan & Stow, 2010). However, 
without direct evidence of greening from ground measurements or 
field observations, conclusions of increased arctic vegetation produc- 
tion from satellite measurements should be viewed with caution. Only 
recently have arctic satellite observations been directly coupled with 
ground observations (e.g. Bhatt et al., 2010), demonstrating that some 
(but not necessarily all) NDVI increases are, in fact, associated with 
enhanced vegetation growth. In our study, direct field observations of 
the site during the study period confirmed that year-to-year changes 
in NDVI during the growing season were attributable to differences in 
productivity and not to differences in cover, standing water, or snow. 
Further field validation of satellite observations will be needed if 
we are to reach unequivocal conclusions regarding arctic vegetation 
change. 

Unlike the ground optical measurements, MODIS was not able to 
clearly resolve interannual differences in snowmelt or NDVI (Fig. 8). 
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While the few “high quality" MODIS NDVI values agreed well with 
field measurements, the majority of the MODIS NDVI values had 
lower quality values, undoubtedly due to the difficulties of “seeing” 
the surface due to the frequent cloud cover following snowmelt at 
this arctic site (Fig. 1A). While it is possible that a more detailed 
analysis of daily MODIS data could improve the MODIS retrievals 
(Narasimhan & Stow, 2010), such an analysis was beyond the scope 
of this study. The ground optical sampling methods reported here 
demonstrate the potential of field optical sampling to provide addi- 
tional insights into the question of the “greening of the arctic” than 
is possible from satellite or gas flux measurements alone. In situ 
field reflectance sampling provides frequent, highly repeatable and 
accurate measurements, avoiding common errors and sampling gaps 
due to clouds. Field measurement allows long sampling periods, en- 
hancing signal-to-noise over satellite sensors, and this is particularly 
relevant under limited arctic illumination. This method can readily 
resolve key events (e.g. snowmelt and rainfall) and microtopographic 
effects on moisture, species composition, and productivity. To better 
clarify patterns of arctic vegetation productivity, we propose systemat- 
ic, automated field sampling of representative arctic locations. Ideally, 
such sampling would be closely coordinated with other measurements 
of productivity (e.g. carbon flux and biomass) and confounding vari- 
ables such as snow cover, standing water, and active layer depth 
would also be monitored to help interpret signals measured from satel- 
lites. The particular advantage of ground sampling is that it can avoid 
artifacts due to cloud cover, which is a frequent summer feature for 
many arctic sites. Given the pronounced effect of microtopography on 
vegetation cover and NDVI, field sampling can also be used to refine 
our understanding of the variables driving vegetation change or stabil- 
ity at a higher spatial resolution than is possible with the coarser 
sampling scale of most satellite measurements. 

5. Conclusions 

Contrary to expectation, earlier growing season did not lead to greater 
greening (NDVI) in this coastal arctic site. These results call into question 
the simple prediction that an earlier snowmelt or longer growing season 
will lead to a greener, more productive arctic. For this site, and other re- 
gions on Alaska's North Slope, recent surface drying trends could cause re- 
duced vegetation production, despite a trend towards earlier snowmelt. 
Clearly, conditions during the growing season, rather than growing sea- 
son length per se, should be considered when evaluating the impact of 
climate change on arctic vegetation production. Unlike field optical mea- 
surements, MODIS satellite NDVI was unable to clearly resolve both the 
timing of snowmelt and the peak season values. We conclude that reports 
of the “greening of the arctic” due to enhanced vegetation indices should 
be taken with caution when based on satellite measurements alone, and 
that ground optical sampling could provide an independent way to eval- 
uate satellite measurements. 
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